XCDCC: Core Excitation in the Breakup of Exotic Nuclei 
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The extended Continuum Discretized Coupled Channel (XCDCC) method is developed to treat 
reactions where core degrees of freedom play a role. The projectile is treated as a multi-configuration 
coupled channels system generated from a valence particle coupled to a deformed core which is 
allowed to excite. The coupled channels initial state breaks up into a coupled channels continuum 
which is discretized into bins, similarly to the original CDCC method. Core collective degrees of 
freedom are also included in the interaction of the core and the target, so that dynamical effects can 
occur during the reaction. We present results for the breakup of "'^'^C^^^C+n and ^^Be^^^Be+n on 
^Be. Results show that the total cross section increases with core deformation. More importantly, 
the relative percentage of the various components of the initial state are modified during the reaction 
process through dynamical effects. This implies that comparing spectroscopic factors from structure 
calculations with experimental cross sections requires more detailed reaction models that go beyond 
the single particle model. 

PACS numbers: 24.10.Eq, 25.60.Gc, 25.70.Do, 27.20.+n 
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I. INTRODUCTION 

Over the last decade, technological developments have 
provided the Radioactive Beam community with detailed 
data on a large variety of reactions involving nuclei far 
from stability Jj. Amongst these are breakup reactions, 
from which properties of loosely bound systems can be 
extracted Q • Breakup reactions have been used to deter- 
mine ground state properties, including binding energies, 
angular momentum structure, momentum distributions 
and spectroscopic factors (e.g. 0, S S S S H ^ 
well as for studying excited state properties and resonant 
states (e.g. [Tol|). 

Breakup is not only of interest for studying the struc- 
ture of these nuclei, but also for astrophysical reasons. 
In a Coulomb dominated process, if we can assume a 
one step semi-classical approximation it is possi- 

ble to factorize the Coulomb dissociation cross section of 
P + t (c + v) + t into a kinematical factor multiply- 
ing the photo-dissociation cross section of the projectile 
c + V P, which in turn is easily translated into the 
capture S-factor at arbitrarily low energies between the 
fragment and core . Manyapplications of this method 
have been used in the pas t fl3 . \lA IT^ and new applica- 
tions are being developed [la |. 

In general, we need to improve on the semiclassical 
models, but still some approximations need to be made, 
depending on the energy regime, the target, and the spe- 
cific kinematical conditions of the experimental setup. 
Recent results show that nuclear contributions extend 
far out into regions where naively one expects they would 
be negligible, because of the extended nature of loosely 
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bound wave functions. Coulomb- nuclear interference is 
also not negligible, the far-field approximation of the 
Coulomb field is inaccurate, and polarization effects need 
to be considered (see for example El El El IM IH ) ■ 
Even when there are no other bound excited states, the 
breakup process goes through several continuum states 
under the influence of the target. These effects have been 
referred to as dynamical effects or continuum-continuum 
couplings, and together with other improvements, have 
helped to solve apparent discrepancies between different 
experiments 23, 24]. 

In earlier days, measurements involving drip line nu- 
clei detected only one of the fragments (e.g. |23). Those 
inclusive data contained inelastic breakup and transfer 
(stripping) as well as elastic breakup (diffraction). In 
the latter, the projectile breaks up and both fragments 
survive after the process without any excitation of the 
target. In stripping, one of the fragments in the projec- 
tile gets absorbed and/or the target gets excited. The 
improved beam intensities in most facilities and the new 
phase of electronics in data acquisition has allowed ef- 
ficient coupling of several types of detectors, such that, 
nowadays, data are less integrated and often contain com- 
plete kinematics (e.g. |2^|23|). For this reason, theoreti- 
cal models need to move toward disentangling the various 
process in their predicted observables. In this work we 
will focus on the elastic breakup component. 

Generally, breakup reactions have been described 
within a three body model, where the projectile is sim- 
plified to a two body system {P = c + v). The relative 
motion between the core and fragment of the projectile 
is distinguished from, and then coupled to, the relative 
motion of the projectile and the target. In the high 
energy regime, straight-line trajectory or adiabatic ap- 
proximations are often made |28ll29l|. At lower energies, 
semi-classical methods taking a Rutherford orbit for the 
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projectile's trajectory are a possible choice fv^ ITM l2^. 
Over all energy regimes, one of the most successful mod- 
els has been the Continuum Discretized Coupled Chan- 
nels method 's^ . It includes the scattering states of the 
projectile completely coupled and is fully quantal. In 
this model, nuclear and Coulomb are treated on equal 
footing, thus interference effects appear naturally. Also, 
any relevant multi-step process within the continuum, or 
from the ground state to the continuum, are automati- 
cally incorporated |^ |^ . 

As successful as it has been, CDCC ^211 has several 
limitations. One of the most serious is the restriction to 
projectiles that can be approximated by two bodies in rel- 
ative motion. The Kyushu group 32] has been expanding 
the original CDCC method to solve problems where the 
projectile is of a clear three body nature, such as Bor- 
romean systems ^^He, "Li, ^^^Be, etc.). Their method 
introduces a Gaussian expansion with a complex argu- 
ment jS.lj ■ essential for computational efficiency. Results 
for nuclear and Coulomb breakup have been recently pre- 
sented [3 1^ • Another limitation is associated with the 
non relativistic description, and this approximation be- 
comes inaccurate at very high energies. Exploratory work 
solving the Klein-Gordon equation [s^ has shown there 
to be non-negligible relativistic effects in the Coulomb 
dissociation of *B for experiments at GSI energies. Fi- 
nally, the projectile-target two-body asymptotics, as an 
approximation of the pure three body asymptotics, has 
also bee n p ointed out as a limiting factor in certain re- 
actions |37| . 

The CDCC method has been applied to many recent 
examples jH [sil S HI [i^. For the ^B breakup 
measured at Notre Dame where only the ^Be fragment 
was detected the theoretical predictions agree very 
well with the data for the kinematic region where only 
elastic breakup is expected jSlj]. If the same ^B struc- 
ture model is used for the breakup of ®B at higher ener- 
gies li^l, one needs to artificially adjust the quadrupole 
strength of the couplings potentials in the reaction [43 . 
This adjustment cannot be understood by a simple renor- 
malization of the E2 strength in the ''Be-p continuum 
p3 | and suggests that structure beyond the single parti- 
cle P3/2 is required to describe the reaction mechanisms 
at higher energies. The elastic and inelastic scattering 
of protons on ^^Be in inverse kinematics was examined 
in the inert core single particle model They found 
that serious discrepancies between the theory and data 
existed in the ^^Be{^°Be(0'^)+n} breakup spectrum in 
the region around resonances built on a 2"*" core. 

There are several ways to improve the two body single 
particle description of the projectile. While the Kyushu 
group has focused on the description of breakup of a three 
body projectile, in this work we will discuss a general- 
ization of the CDCC method for a two body projectile 
where the wave function is no longer single particle. This 
generalization of the CDCC method, which we denomi- 
nate extended CDCC (XCDCC) 46], includes a cou- 
pled channels description of the projectile's initial and 



final states, but also any dynamical core excitation/de- 
cxcitation that may occur in the reaction process. 

It is clear that if core excitation/de-excitation happens 
during the reaction process, this will have implications for 
not only the excited cross sections, but also for measur- 
ing the ground state core distributions. In addition, one 
expects that the excitation/de-excitation process will de- 
pend strongly on the type of reaction, namely high versus 
low energy, light versus heavy target. Even when there is 
only a single valence neutron/proton, the valence nucleon 
may occupy several configurations, in some of which the 
core is excited. By measuring the gamma rays in coinci- 
dence with outgoing fragments, core states can be iden- 
tified, and spectroscopic factors for the overlap of the 
projectile ground state with specific core states can be 
examined |47l ]. 

In nuclear dominated knockout reactions on light tar- 
gets, typically the method of ^3] is followed. Therein 
one calculates the single particle cross section for both 
diffraction (elastic breakup) and stripping, in an eikonal 
model, for each ground state configuration of the projec- 
tile. The single particle cross sections for each state are 
then multiplied by the spectroscopic factor for that con- 
figuration, calculated within a shell model, and the total 
knockout cross sections is obtained by adding the sin- 
gle particle contributions incoherently. Improvements on 
the eikonal approximation for the breakup contribution 
E^li^ have been examined using CDCC. These models, 
however, are still single particle and neglect any inter- 
ference between configurations, as well as any dynamical 
excitation of the core in the reaction. We will examine 
the importance of these effects in this paper. 

In Coulomb dissociation, by contrast, a more simplistic 
approach is typically assumed: that the relative probabil- 
ities do not change during the reaction, and the contribu- 
tion from excited core components is negligible 9] . Thus 
the ratio of the theoretical cross section, for the ground- 
state-core configuration, to the experimental cross sec- 
tion is often taken to be directly the spectroscopic factor 
for the ground state. 

Contributions from excited core components in 
Coulomb dissociation have previously been estimated 
|5lj ]. by calculating the contribution from excited core 
components separately, and adding them incoherently, 
following the method of ^^l- The Coulomb dissocia- 
tion was calculated using DWBA and Adiabatic reac- 
tion models. In [HJ], it was concluded that the Coulomb 
dissociation to excited core states was very small, only 
a few percent, very much less than in nuclear breakup. 
They then suggest that since the cross section from ex- 
cited core states has a negligible contribution to the total 
Coulomb breakup cross section, one could compare the 
single particle theoretical cross section to the experimen- 
tal cross section to obtain the spectroscopic factor for the 
core in its ground state. One has to consider the many 
approximations assumed in the theory. This model ne- 
glected any interference between the configurations and 
any dynamical (de-)excitation of the core. These effects 
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can be examined with XCDCC, but Coulomb dominated 
reactions will be a subject of a future publication. 

One simple way to introduce collective degrees of free- 
dom to the core is to assume it behaves as a perfect ro- 
tor. Rotor-\-nucleon coupled channels models for loosely 
bound systems have been a ppli ed to a number of light 
nuclei on the dripline (e.g. |52L IsM \5^). This type of 
structure model for the projectile has been also used in 
the calculations of reaction observables [H^ [s^ . 

In [s^ , the coupled wave functions are used to compute 
the B(E1) electric excitation function from the ground 
state into the continuum. The interference between the 
states with the same core spin in a pure Coulomb disso- 
ciation process is then studied within the semi-classical 
first order approximation JsS^. Interference effects be- 
tween states with different core spins were neglected since 
this model does not include inelastic transitions caused 
by the target. These interference effects were shown to 
be important in ^^C on 208pi3 gy ^ MeV. 

In 56] , the eikonal approximation is used to calculate 
the total cross sections for the breakup of a deformed core 
plus nucleon system. Results for ^^C breakup on ^Be at 
62 A MeV show a dramatic increase of the cross section 
with the quadrupole deformation parameter. Contrary 
to the semi-classical first order approximation, this model 
includes dynamical effects, but only total cross sections 
are obtained. It can be interpreted as a precursor to the 
new model presented here. 

There is a wide variety of problems that can be inves- 
tigated with XCDCC, and this is the first chapter where 
the formalism is given in detail and integrated observ- 
ables for two cases are presented, namely ^^C and ^^Be 
breakup on ^Be at « 60 A MeV, both measured at MSU 
0, • We present a comparison with the data of ?| in a 
shorter paper [i^ . 

In section the theoretical development is presented. 
A description of the structure models used for the pro- 
jectiles is detailed in section IIIII Numerical details for 
the examples considered are given in section ITVI Follow- 
ing, in section the results are detailed and discussed. 
Finally, section IVll contains the conclusions drawn from 
our study and an outlook for future work. 

II. THEORETICAL DEVELOPMENTS 

The breakup of a projectile P consisting of a valence 
particle (v) loosely bound to a core (c) on a target t can 
be modelled as a three body scattering problem. The 
three body Hamiltonian for the system has the form: 

H = Tr + i/p,oj + Vet + Kt , (1) 

where represents the kinetic energy operator for the 
projectile-target system, iJproj is the internal Hamilto- 
nian of the projectile and Uct {Uvt) is the interaction 
between the core (valence) particle and the target. The 
Hamiltonian is usually written in the Jacobi coordinates 
(i?, r) where R is the coordinate for the center of mass of 




FIG. 1: Coordinates for the three body scattering problem 
to be considered. 

the projectile relative to the center of mass of the target, 
and r is the coordinate of the valence particle relative to 
the center of mass of the core as shown in Fig. ^ 

The full three body wave function for the system with 
definite total spin Jt is expanded over the states of the 
projectile, both bound and unbound: 

= E <^ (^) [[^i (^) ® ^)]^ ® (^*)] t ■ 

a. 

(2) 

The XCDCC states, labelled by a = {L, Jp, J, Jt,i,n}, 
describe the projectile state {Jp,i,n}, where Jp is the 
total spin of the projectile, i is the discretization in rel- 
ative momentum and n is an additional boundary con- 
dition that is required to describe coupled channels bins 
and is discussed later. In this paper we assume the tar- 
get is structureless, therefore removing the dependence 
on (^t) in Eq. The projectile-target orbital angular 
momentum, L, and the spin of the target Jt all couple 
to a total spin of the three-body system Jt- We use the 
notation 

(fi,r,^|a; Jt) = [[¥^(11) ® <PJAr, Oh ® <i>J.] (3) 

to describe the angular momentum and projectile wave 
function compactly. 

If the interaction between the fragments in the projec- 
tile depends only on r, the projectile Hamiltonian -ffproj 
generates single particle bound state wave functions, and 
single particle continuum states which can later be dis- 
cretized into continuum bins. In this work, we wish to 
include an additional degree of freedom, ^, concerning 
the projectile's core, and allow the core states to be cou- 
pled together. For simplicity we take a rotational model 
for the core |53 | , but the evaluation of the coupled chan- 
nels bins and potential matrix elements are not specific 
to the choice of this collective model. We consider that 
the core-valence interaction depends on core collective 
degrees of freedom ^: 

Hproi =T, + Kc(r, + /jcorc(C) , (4) 

where T^. is the core- valence kinetic energy operator, and 
hcoreiO is the intrinsic Hamiltonian for the core. We de- 
fine (fi to be the eigen function of the core corresponding 
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to eigen-energies ej. Now the projectile ground state con- 
tains contributions from different core states, in a fully 
coupled manner: 



'>Xs\ 



(5) 



We use the subscript on the projectile wave function to 
denote the ground state (the supercripts (i, n) in Eq. (jS)) 
are superfluous for bound states). Here, a given chan- 
nel a depends on a set of quantum numbers describing 
the relative (/) and total (j) angular momentum of the 
valence particle, with fixed spin s, which couples to the 
core spin I. We label each channel which can couple to 
total spin Jp by a = {l,j,I}. We use the notation 



Jp 



(6) 



to describe the angular momentum and core wave func- 
tion compactly. Previous work ^52j gives details of such 
rotor +n models. 

In addition to the dependence of the internal Hamil- 
tonian of the projectile on ^, the core-target interaction 
Vct{R,£,) also contains the dependence on the degrees of 
freedom of the core and therefore, for a breakup pro- 
cess within this model, the core-target interaction can 
excite/de-excite the core during the reaction process. 
The coupling matrix elements of this deformed potential 
will be discussed in section Hi Bl 



A. Construction of projectile coupled channels bins 

The projectile's {c-\-v) scattering states for a given in- 
coming kinetic energy Ek are solutions of a coupled chan- 
nels equation which can be reduced to the following radial 
form: 

[Ek-Sa+en-Ta]fa:n{r;kan) = '^Va:a' (r) fa' -.nir] ka' n) , 

a' 

(7) 

where a and n stand for outgoing and incoming bound- 
ary conditions respectively, — /i^fc^/2/i^c is the core- 
valence relative energy, k is the core- valence relative mo- 
mentum, and the reduced mass is fj,vc- The momentum 
that appears in the scattering wave, fcan; Is the asymp- 
totic momentum available in particular channel an, and 
is defined by 



ft2 



(8) 



where Eq is the energy of the core in the outgoing chan- 
nel a, and e„ is the energy of the core in the incoming 
channel, n. The kinetic energy operator, 

is defined for the angular momentum I of the core- valence 
system in channel a. The coupling potential matrix is 

K:a'(r) = {{ls)jJ;Jp\VUr,0\{l's)j',I';Jp) , (10) 



(9) 



where a = {l,j,I} are the outgoing channel quantum 
numbers. 

The open scattering wave functions arc obtained by 
matching to the asymptotic boundary condition. 



fa:n{r; k) = -[H-{kr)5an " i?+(fcr)5,^-l 



(11) 



at rmatch to obtain the coupled channels square S'-matrix, 
Sir^. Here Hi^\kr) are the diagonal Coulomb func- 
tions. Inelastic channels with excited core components 
will remain closed unless the relative energy is higher 
than the core energy or the incoming channel has an ex- 
cited core component. If the available energy for channel 
(an) is negative, then the wave function fa.nir; kan) will 
be closed with separation energy ?i^/c^„/2^^c- 

In principle, coupled channels scattering wave func- 
tions could be used directly to calculate the one-step 
breakup process. However, as multi-step effects can be 
extremely important, we opt to discretize the coupled 
channels set, into square integrable bin wave functions. 
The generalization of the standard binning procedure js^ 
to coupled channels bins is described below. 

For each projectile Jp state in the continuum, we de- 
fine coupled channels bin wave functions that are square 
integrable, resulting from a superposition of pure scatter- 
ing states (calculated from Eq. Q) within a momentum 
bin [ki-i,ki\. Each momentum interval is labelled by i. 
The coupled channels bin wave function are written as 



ttNI 



ki 



dk gn{k)fa:n{r;kan) (12) 



where, as before, k is the relative momentum between 
core and valence fragments. Here, fco represents the 
breakup threshold between the core and valence parti- 
cles. 

The weight factor of the coupled channels bins, g-nik), 
is chosen to be e""'"'-'^-', where (S„(fc) is the phase shift of 
channel a = n (i.e. the phase obtained from the diagonal 
of the S'-matrix). In the one channel case, this phase 
makes the wave functions real '2l'| . The coupled channel 
bins are now complex, and the full complex wave function 
is used in the CDCC equations up to the matching radius 
''match- By choosing the weight factor as the phase of 
channel a = n, the wave functions are approximately real. 
This is useful when matching to the asymptotic Coulomb 
functions '5^ in the CDCC equations, as the long-range 
Coulomb functions become approximately real. 

The normalization factor is 



\9n{k)\'dk 



(13) 



Since our weighting has |g„(A:)p = 1, the normalization 
is simply = ki — fci_i, and the avera ge e nergy of bin 
a is = h^{kl - kf_;^)/6fi.vc{ki - fc^_i) llj. 

The full bin wave function describing the projectile 
state Jp, for momentum bin i and incoming wave n, 
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is the sum of the coupled channels bins obtained from 
Eq. 



(14) 



The incoming channels form an orthogonal set and there- 
fore projectile wave functions carry the label n. It is im- 
portant to note that since the bin wave functions appear 
in the bra in the three-body T-matrix, these incoming 
wave boundary conditions on the scattering wave relate 
to the outgoing channel of the breakup fragments. 



[53, by the theory of Appendix IbI The new physical ele- 
ment is in the core-target interaction, which now depends 
on the core degrees of freedom (^). 

The core-target interaction can be expanded in multi- 
poles of Q, 

V,t{R.r,0 =Y.^^^S{r.)Y.^Q,{r<^)C*Q,{S,) (18) 
Q 



q 



where = R—jr is the core-target coordinate, with 7 = 
niy/mp as the ratio of the valence mass to the projectile 
mass, and CQq{fc) is related to the spherical harmonic. 



B. The three body coupling potentials 

The scattering problem for the three body Hamilto- 
nian, Eq. ^ H'^ = iSvP, can be reduced to a set of 
radial coupled channels equations in R when using the 
expansion Eq. (O . Then the scattering coupled channels 
equation becomes 

(15) 

where the kinetic energy operator is 



Tr 



2fi 



pt 



L{L + l) 



di?2 



i?2 



(16) 



and /ipt is the reduced mass of the projectile-target sys- 
tem. The energy of each channel = -E — is the 
total energy of the three-body system minus the aver- 
age energy of bin a. The incident center-of-mass energy 
will he Eq, as a = Q represents the ground state of the 
projectile. 

The coupling potentials for the elastic breakup of the 
projectile (the process where the target remains in its 
ground state) are 



UiT^,{R) = {a-:JT\V,t{R.r,0 



VMr)\a';JT) , 

(17) 



where a describes the angular momentum as defined in 
Eq. The projectile-target interaction potential is the 
sum of the core-target [Vet] and valence-target [Vyt) po- 
tentials, and contains both nuclear and Coulomb terms. 
The valence- target interaction can be treated as before 



(19) 



where we use the notation Q = ^/2Q + 1. 

Next, the product r^CQq{rc) can be expanded in mul- 
tipoles A for the r and R coordinates, as in [s^: 



r?CQq{r,) =Q^E 



(20)! 



(2A)![2(g-A)] 



■i?^(-7r)' 



C\^{R)Cq- 



\,q-fj. 



(r) 



-A A g 



(20) 



The multipole Q for the core-target potential (re) 
can be expanded in multipoles of 



K?(rc) 



Y,k^V^,''{r,R)Y,CKu{R)Cl^{r) , (21) 



'■^ K k 

where the radial potential dependence is given by 



— FKWau ; u = Hr. 

L 

(22) 

Combining Eqs. 10, HOI) and we find that the 

spherical harmonics in R can be coupled to a new multi- 
pole A and, similarly, the spherical harmonics in r can be 
coupled to another multipole A'. The angular momen- 
tum structure gives three spherical harmonics in R, r, ^ 
coupling to form a tripolar spherical harmonic of rank 
zero. Summing over the angular momentum projections 
(q,/x. A;), the core-target potential operator can then be 
reduced (according to Appendix to 



V.tiR,r,0 = E vS''{r,R)R\~ir)' 

KQX 

X E AA'' 

AA' 



(20)! 



(2A)![2(0-A)] 





A 




(0 












A 0-A K 



(23) 



{CA{R)-{CA,{r)<E>CQ{0}A) ■ 
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The transition couphngs of Eq. H17() are calculated as 
matrix elements of Eq. (|23l) , assuming the coupling order 
defined in Eq. It is useful to write the transition 

couplings in terms of a transition potential which only 
depends on the overall projectile initial and final states, 
J'p, i' , n' and Jp, i, n, and the transition multipolarity A. 
Combining the potential operator from Eq. (|23|l and the 
matrix elements from Appendix ^ we write the core- 
target three-body couplings as 

U'JAR) ^LL'Jpjli-iy-+-'J2i-l)^A^ 





L 




(Jp 


J'p 













L 





J^Jpin-.J'pi' n' i^) • 

(24) 

The form factor is the sum over KQX multipoles and the 



projectile coupled channels (a and a'), 

KQX 

where the radial dependence is defined by the integral: 

r,KQ\ / r}\ 

^ain:a'i'n' V^/ 

= K r"\l^{rrvS''{r,R)R\-^r)Q-\i,.,^,{r)dr . 
Jo 

(26) 

The remaining angular momentum dependence is sum- 
marized as 



^z^^l^o ojl^oooJ\Ag-Ai^/V'/s/1jj, ^^r 



The angular momentum coupling in Eq. H27I) contains 
the matrix element (/ || Cq{£) \\ I'), which describes 
elastic and inelastic excitations of the core through its 
scattering by the target. In this paper we use a rotational 
model for the core, for which this matrix element has a 
simple solution (Appendix Q. 

This matrix element (/ || Cq(^) || /') is left separate in 
Eq. (|27|l since, in general, other models could be used to 
induce the inelastic scattering of the core. Given a core- 
target interaction with core degrees of freedom, ^, and 
multipole expansion, Q, the inelastic modes can be calcu- 
lated, and along with radial coupling potentials [Eq. (|10|) ] 
for the core- valence interaction, XCDCC opens the door 
for other models of the core which go beyond the collec- 
tive model. 

The monopole couplings for the core-target interaction 
[Q = in Eqs. (|24I27|) ] are evaluated in Appendix |b1 We 
see by comparing Eq. (|24|l and Eq. HB1|I that including 
extra core degrees of freedom for the core does not in- 
crease the complexity of the final form factors, since they 
both depend only on a multipolarity and the transition 
J'p — > Jp. Thus the form factors can be pre-calculated 
efficiently. To reduce memory requirements, the individ- 
ual components Rain^i'i'n' iR)PaS^^ '^^'^ ^c summcd so 
only the final Ffpi^-j' i'n'(^) stored. A new version of 
Fresco jsl] was created for this purpose. 

The only increase in the size of the calculation com- 
pared with the single particle calculations, with no core 



degrees of freedom, is the increase in the model space 
due to the extra quantum number n for the projectile 
state. This relates to the incoming wave boundary condi- 
tion on the coupled channels solution of the core- valence 
wave function. One needs to include all possible incoming 
waves for each of the coupled channels scattering states. 
This extra quantum number carries the information re- 
garding the final state of the core-valence system in the 
CDCC equations. 



III. STRUCTURE OF ^^BE AND ^^C IN THE 
ROTATIONAL MODEL 

We use a rotational model to describe the projectile 
as a coupled core- valence system. The rotational model 
for ^''C is based on the model used in [s^l- We use the 
same potential parameters as 0, given in Tabled We 
neglect the spin of the valence particle, because we are 
interested in the effect of the core spin, and the effect 
of the valence particle spin on the total cross sections 
is small. This model produces a ^''C ground state spin 
of Jp ~ 2+. We use a quadrupole mass deformation 
of (32 = 0.55, as in [5^. Modeling the ^^C core with two 
states {0"*", 2+}, with an excitation energy of E2+ = 1.766 
MeV, and including all possible relative angular momen- 
tum up to ^ = 4, produces four coupled channels in the 
ground state: [d(g)0+], [s(g)2+], [d®2+], [g®2+],wii\\ 
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model 


interaction 


V Rv av /32 


CC 
SP 
SP 
SP 


^''C{0+,2+}+n 
^«C{2+}+n{d} 
i«C{2+}+n{s} 
^«C{0+}+n{d} 


62.0 3.2 0.65 0.55 
52.7 3.2 0.7 - 

49.4 3.2 0.7 - 

48.5 3.2 0.7 - 



TABLE I: ^''C+n potential parameters for Eq. dHTll . 

relative probabilities of 0.175, 0.103, 0.721, 0.001 respec- 
tively. This model predicts the dominant channel to be 
[d® 2+], which is consistent with the large spectroscopic 
factor obtained from shell model calculations (this chan- 
nel was left out of the model space in the rotational model 
of 13). The single particle (SP) central potential pa- 
rameters are taken from |56j |. and are adjusted to fit the 
ground state binding energy for each configuration. ^^C 
also has other bound states around 4 MeV in energy, one 
of which, the 4+ at i?4+ = 4.142 MeV, could be included 
within our rotational model. In some of the calculations, 
we include this core state in the projectile continuum al- 
though we keep the projectile's ground state description 
without the 4+ core state for consistency. Using a rms 
matter radius of 2.70 fm [g^I for the ^^C core, the rms 
radius calculated for the deformed ^^C is 2.77 fm. The 
rms radius for the single particle [d (S) 2+] state is 2.78 
fm. The weighted sum for the incoherently summed sin- 
gle particle projectile states, has a rms radius of 2.79 fm. 
The experimental determined rms radus is 3.04±0.11 fm 

For ^^Be we base our potential parameters on those of 
[s^ . These were fitted to give the correct binding ener- 
gies of the ground and first excited state, as well as the 
B(E1) transition between them, using a parity dependent 
potential. Since again we neglect the spin of the neutron, 
and therefore no spin-orbit force, we readjust the poten- 
tial parameters slightly to fit the binding energies of the 
bound states. The parameters used in the calculation 
are given in Table (nj with the model CC referring to the 
coupled rotational model for the projectile, where a de- 
formation P2 = 0.67 is taken from [53 ■ This model pro- 
duces a ^^Be ground state with spin Jp = 0"*", with two 
coupled channels: a s-wave neutron coupled to O"*" core, 
and a d-wave neutron coupled to a 2+ core. The rela- 
tive probabilities are 0.883 and 0.117 for the [s(8)0"'"] and 
[d(g)2"'"] states respectively. The single particle model for 
the core consists of a parity dependent potential as used 
in [s^ , which fits the bound states using a central poten- 
tial. For the states with a 2+ core the potential depth is 
adjusted to fit the binding energy plus the energy of the 
2+ core. Using a rms matter radius of 2.28 fm [g^ for 
the -'^"Be core, the rms radius calculated for the deformed 
"'^"Be-l-n is 2.99 fm. The rms radius for the single particle 
O"*" core calculation is 3.08 fm. The weighted sum for the 
incoherently summed single particle projectile states has 
a radius of 3.01 fm. These compare reasonably well to 
the rms radius for "Be of 2.90±0.05 fm |6||. The radii in 
our few-body model are a bit small, but since we do not 
compare to data, and only compare to previous work, we 



model 


interaction 


parity 


V 


Rv 


av 


132 


CC 


^"Be{0+,2+}+n 


-1- 


55.25 


2.483 


0.65 


0.67 


CC 


"Be{0^, 2+}+n 




47.00 


2.483 


0.65 


0.67 


SP 


i^BejO+j+n 


+ 


55.50 


2.736 


0.67 




SP 


'"BejO+l+n 




30.48 


2.736 


0.67 




SP 


i''Be{2+}+n 


+ 


75.07 


2.736 


0.67 




SP 


iOBe{2+}+n 




39.95 


2.736 


0.67 





TABLE II: ^"Be+n potential parameters for Eq. dHil . 



keep the same potential model for consistency. 



IV. NUMERICAL CONSIDERATIONS 

In all the calculations the core-valence integrals were 
performed out to a maximum radius of rmatch = 70 fm, 
with a step length oi h = 0.05 fm. The momentum in- 
tegrations for each bin [Eq. 1)12(1 ] were calculated with 
60 steps. The multipole transitions were calculated with 
all possible K, Q and A transitions, with the transition 
multipolarity truncated to A < 4. Partial waves were 
included up to L = 5000 for ^^C breakup and up to 
L — 3000 for ^^Be breakup. The radial wave functions 
were matched to the asymptotic Coulomb functions at a 
radius of i?asym = 500 fm. The spin of the valence par- 
ticle was neglected and therefore the channel labels all 
have j = I. The spin of the ^Be target was set to zero 
as it is a spectator in our reaction model. The channel 
labels refer to the total spin and parity of the projectile 
(Jp), the angular momentum of the neutron (Z), and the 
spin of the core fragment (/), in the final state. 

The XCDCC model space for the ^^C reaction where 
i^C is modeled with two states {0+,2+}, the XCDCC 
model space is given in Fig. |21 For the -'^^C reaction 
where ^^C is modeled with three states {0+, 2+, 4+}, the 
XCDCC model space is given in Fig. E| The XCDCC 
model space for -'^^Be breakup is given in Fig.^ The num- 
bers of bins spanning an energy range for each channel in 
the final state are given in each box, with the bins at reg- 
ular intervals in momentum. Summaries of the various 
sizes for the model space are given in table IIIll The first 
column is the projectile model space: which core states 
are included in the calculation. The next five columns 
give the number of permutations for each of the variables; 
the number of coupled channel states {a J pin), the num- 
ber of projectile states {Jpin), the number of couplings 
between coupled channel states {P^^^'^), the number 
of form factors stored {Fj^-^.j, j,„,), and the number of 

channels in the CDCC equations (a). The deformed core 
calculations were performed on a SGI Altix 3700, while 
the single particle calculations were performed on a clus- 
ter of AMD Opteron 250 processors. 
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Projectile 



aJpin Jpin P 



KQ\:A 



CPUtime(s) NCPUS walltime(h:m) Mem(GB) 





90 




7312 


258 


20568 


11 


00:35 


0.2 


572 


179 


3862133 


58951 


542 


342008 


8 


13:32 


1.2 




103 




9788 


344 


37900 


1 


10:13 


0.2 


478 


144 


2886779 


41795 


473 


206048 


15 


03:59 


1.4 


1269 


244 


40036229 


118450 


546 


2991231 


16 


56:01 


2.5 



^"Be{0+}+n 
i°Be{0+,2+}+n 

i'^C{2+}+n 
i'5C{0+,2+}+n 
^^C{0+,2+,4+}+n 



TABLE III: Size of calculations for the different model spaces. The first column gives the projectile model, with the core states 
included. The next 5 columns give the number of permutations for various variables that arc required. The last 4 columns give 
total CPU time, the number of CPUs, walltime, and memory requirements to store the form factors. 



type 


interaction 


V rv av W rw aw 


central 
deformed 

central 
central 


core + target 
core + target 

neutron + target (E=62 MeV) 
neutron + target (E=60 MeV) 


123.0 0.75 0.8 65.0 0.78 0.8 

134.1 0.75 0.8 68.15 0.75 0.8 
34.54 1.17 0.75 13.4 1.26 0.58 
33.79 1.17 0.75 12.08 1.26 0.58 



TABLE IV: Fragment-Target optical potential parameters. The radius parameters are to be multiplied by AcoL + ^target fc"" 
the core-target potential and ^target fo'' the neutron-target potential. The core-target parameters are used for both ^^C and 
'^"Be, and the neutron-target parameters are only adjusted slightly for the different energy of the reaction. 
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FIG. 2: ^'^C continuum model space. The number of bins and the energy range are given for each outgoing channel {1,1^) for 
each spin parity combination of the projectile (Jp). 
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FIG. 3: '^C continuum model space. The number of bins and the energy range are given for each outgoing cliarmcl {I, I^) for 
each spin parity combination of the projectile (Jp). 
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FIG. 4: ^^Be continuum model space. The number of bins and the energy range are given for each outgoing channel I'^) for 
each spin parity combination of the projectile (Jp). 



V. RESULTS 

Here we present the results of the calculations. For 
both reactions, we perform three different calculations 
to understand the relevance of the different ingredients. 
First we simulate the CDCC version of the single parti- 
cle calculations in [^Q- We refer to it as SPIS: single 
particle incoherent sum. It corresponds to a sum over 
the CDCC single particle cross sections multiplied by the 
corresponding relative probability for the projectile gen- 
erated in Sect ion UTTl It does not contain any interference 
between core states nor core excitation/de-excitation due 
to the interaction with the target. Secondly, we perform 
a calculation where we do not include dynamic effects 
with the target, and refer to it as CCSE: coupled chan- 
nels static excitation. It contains the full coupling within 
the projectile but has only Q = in the three-body cou- 
pling matrix elements. The third case is the full calcu- 
lation, which we denote by CCDE for coupled channels 
dynamic excitation. In all the above, we perform the 
calculations switching on and off continuum-continuum 
coupling to determine the importance of these couplings 
as done before In Tables IVHVIIII the couplings la- 
belled c-v and c-T refer to the deformed couplings be- 
tween the core-valence and core-target subsystems. The 
monopole couplings are always included. 

A. ^Be(i''C,i^C+n) 62 MeV/nucleon 

The 1-step (no continuum-continuum coupling) 
breakup cross sections for ^Be(^^C,^^C{0+, 2+}+n) are 
given in Table IVI The single particle cross sections for 
each of the configurations used in the CC calculations 
are 33 mb, 25 mb and 13 mb for the [d (g) 0+], [s (g) 2+], 
and [d (g) 2+] states respectively (we neglect the [g (g 2+] 
state since its contribution is very small) . These are then 
multiplied by 0.175, 0.103, and 0.721 to obtain the single 
particle incoherent sum SPIS cross section of 18 mb. This 



Couplings 

Model c-v c-T (Tg+ (T2+ <j 

SPIS '~ X 6 12~ T8 

CCSE / X 5 11 17 

CCDE / / 7 14 21 



TABLE V: "C 1-step (no CC couplings) breakup cross sec- 
tions leading to ^''C{0+,2+} + n: SPIS is the single particle 
incoherent sum model, where neither core-valence or core- 
target excitation couplings are included, CCSE is the coupled 
channel static excitation calculations, where core-valence cou- 
plings are included but no core-target, and finally the full cal- 
culation CCDE coupled channel dynamic excitation, where 
both core-valence or core-target excitation couplings are in- 
cluded. 



model is equivalent to that used widely to analyze knock- 
out data to specific states difference being that 
this is a 1-step breakup reaction model which includes nu- 
clear and Coulomb breakup whereas the eikonal model 
includes single-particle-continuum-continuum couplings 
but nuclear breakup only. We will include continuum- 
continuum couplings in our model later. Introducing a 
coupled channels description of the projectile (CCSE) 
provides very similar cross sections when comparing to 
the SPIS model. This suggests little interference between 
configurations during the reaction. Including dynami- 
cal excitation in the reaction model (CCDE) enhances 
the cross section to both core states, increasing the to- 
tal breakup cross section to 21 mb in agreement with 
the findings of Batham et al. 56j. However in that 
work, only cross sections summed over all core states were 
given. 

We explore the relevance of other states in ^^C by first 
calculating the inelastic reaction ^Be{^^C{2+),^^C*). We 
use a simple rotational model for ^^C, and include three 
states {0+, 2+, 4+}. We have ^^C{2+) in the initial state 
since "'^^C has predominantly ^^C(2+) in the ground state. 
The resulting inelastic cross section to specific core states 
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Couplings 










Model 


c-v 


c-T 


CC 




(72+ 




O" 


SPIS 


X 


X 


X 


6 


12 




18 


CCSE 


/ 


X 


X 


6 


12 


5 


22 


CCDE 


/ 


/ 


X 


7 


15 


5 


27 


SPIS 


X 


X 


/ 


5 


11 




16 


CCDE 


/ 


/ 


/ 


4 


10 


5 


19 



TABLE VI: ^^C breakup cross sections leading to 
^'^C{0+, 2+, 4+} + n. The abreviations in the first row are the 
same as in Table IVI Calculations are repeated with and with- 
out continuum-continuum couplings (labeled here by CC). 





Couplings 








Model 


c-v 


c-T 


CC 




0"2+ 


a 


SPIS 


X 


X 


X 


140 


1 


141 


CCSE 


/ 


X 


X 


128 


7 


135 


CCDE 


/ 


/ 


X 


142 


8 


150 


SPIS 


X 


X 


/ 


109 


1 


110 


CCSE 


/ 


X 


/ 


107 


8 


115 


CCDE 


/ 


/ 


/ 


109 


8 


117 



TABLE VII: ^^Be nuclear and Coulomb breakup cross sec- 
tions for ^°Be{0''', 2"*"} + n The abreviations in the first row 
are the same as in Table IVI 



is CTo+ = 4 mb, CT2+ = 5 mb, and (T4+ = 12 mb. The 
large number of projections of the 4+ state and its rela- 
tively low energy = 4.142 Mev, raises suspicion that 
this 4+ state will hold a significant contribution to the 
breakup cross section of "'^^C. We recognize that the 2+ 
excited state in -'^^C is not described well in a collective 
model, but here we use it as a toy model to compare 
with previous work. One may also wonder about contri- 
butions from other states in a rotational band, since we 
are using the particle-rotor model. We include a 4+ state 
here since we start with ^^C in a 2+ state. Inelastic exci- 
tations with A/ = 4, within our rotational model, are an 
order of magnitude smaller than the A/ = 2 transitions, 
thus it should only be necessary to include the A/ = 2 
states for the core in our rotational model. 

We therefore recalculate the cross section for ^^C 
breakup including the 4+ state in the continuum (Ta- 
ble IVip . We keep the same ground state with only 
{0+, 2+} so that the relative weights of each ground state 
configuration are preserved. Since we have no 4"*" in the 
ground state, the SPIS model is the same as before. The 
4+ state enhances the total breakup cross section in the 
CCSE model from 17 mb to 22 mb and in the CCDE 
model from 21 to 27 mb. The last two rows in Ta- 
ble IVI! are calculations including continuum-continuum 
couplings. These couplings suppress the cross section as 
was seen before in the single particle case ^|. It 
is a mere coincidence, though, that the final total cross 
section of 19 mb compares so well with the SPIS result. 



B. ''Be("Be,"Be+n) 60 MeV/nucleon 

The breakup cross sections for the reaction 
9Be("Be,i°Be{Q+,2+}+n) are shown in Table IVni 
The single particle cross sections for 1-step breakup are 
159 mb for the [s ® 0+] state, and 10 mb for the [d ® 2+] 
state. Multiplying these single particle cross sections by 
the relative probabilities of 0.883 and 0.117 respectively, 
one gets a total breakup cross section of 141 mb for the 
SPIS model. If we calculate the cross section in the 
CCSE, the cross section drops to 135 mb, suggesting 
some destructive interference between the configura- 
tions. When we include dynamical excitation (CCDE) 
the cross section is enhanced at 150 mb. The apparent 



agreement between the cross section to the 0''^ state (142 
mb) in the CCDE model and the single particle cross 
section multiplied by the relative probability (140 mb) 
is merely a coincidence. It is a result of cancellations 
between two different effects, the interference between 
different configurations of ^-'^Be and the dynamical 
excitation couplings in the core-target potential. We 
also see that the SPIS model severely underpredicts the 
amount of cross section to the 2+ state. The additional 
cross section seems to come from including the couplings 
between the ^°Be and the neutron rather than from 
the dynamical excitation couplings of the core-target 
interaction. When continuum-continuum couplings 
are included there is a general reduction of the cross 
sections, but the same qualitative effects of projectile 
couplings and dynamical excitation remain. 

In the eikonal model of breakup reactions used in the 
analysis of knockout data [4^ . only nuclear breakup is 
included. In order to compare with previous calcula- 
tions for this reaction 0, where the contributions from 
Coulomb breakup and dynamical excitation were roughly 
estimated, we rerun our calculations with only nuclear 
breakup. Result for the CCDE show that the nuclear 
cross section to the 0+ state is 10 mb less than the cross 
section calculated when both nuclear and Coulomb are 
included. This appears to agree well with the estimate 
of 10 mb given in [3], however it does not mean that 
Coulomb and nuclear should be added incoherently. A 
Coulomb breakup only calculation yields 12 mb. The 
estimated contribution from dynamical excitation of 0"*" 
core to the 2+ excited state was 11 mb multiplied by the 
spectroscopic factor of the 0+ core, equivalent to 0.883 
in our rotational model, which gives 10 mb. We see from 
Table IVIII that the increase in the cross section to the 
2+ state, from the SPIS to our CCDE model is 7 mb. 
So it turns out that the estimates presented in ^7] pro- 
duced total cross sections close to those we obtain in a 
more accurate calculation. Again, it is important to note 
that in general this will not happen. A comparison of 
the knockout data from within this new model which 
incorporates core deformation in a consistent manner is 
given in |4^ . 
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/ 
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7 
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TABLE VIII 
^°Be{0+,2+} + n 
same as in Table IVI 



^'^Be nuclear breakup cross sections for 
The abreviations in the first row are the 



VI. 



SUMMARY AND CONCLUSIONS 



In summary, we have extended the CDCC method to 
inchidc a coupled channels description of the projectile, 
using a deformed core-valence interaction, thus allowing 
excited core contributions in the projectile states. We 
also include a deformed core-target interaction, allowing 
the core to (de-)excite during the reaction. 

For the reactions considered here, the breakup of ^^C 
and ^^Be on a ^Be target at « 60 MeV/nucleon, the 
comparison of the structure calculations with the experi- 
mental data will depend on the reaction model used. By 
including more of the important dynamics of the reac- 
tion mechanisms, a more direct comparison of structure 
and experiment will be possible. Even though the ex- 
cited states of ^^C are not well described by collective 
excitations, we use a toy model of rotational excitations 
to the excited 2+ and 4"*" states. This is to compare with 
previous calculations. In the breakup of "'^^C we find that 
total cross section is significantly increased due to ex- 
citation of the core, which is consistent with previous 
calculations. Our calculations here go beyond previous 
work allowing for a partial wave decomposition of the 
final core states. Even when there is no 4+ core-state 
in the ground state of ^^C, significant population of the 
4"*" is predicted in the final breakup states. Continuum- 
continuum couplings decrease the cross section. 

The particle-rotor model for ^^Be describes the nucleus 
well, and there is good agreement with the spectroscopic 
factors obtained from shell model calculations. When we 
allow the ^"Be to deform, we see an increase in the total 
breakup cross section compared to the weighted incoher- 
ent sum of the single particle cross sections. The cross 
section to the 2+ state in ^°Be is significantly increased 
when the core is deformed. 

From our calculations we conclude that both the de- 
formed couplings between the core and valence particles 
and the deformed couplings between the core and tar- 
get are important and should be included in the reaction 
model. This eXtended CDCC (XCDCC) reaction model 
can now be used to describe breakup reactions where sig- 
nificant contributions from excited core states are needed 
to describe the projectile at any point during the reac- 
tion process. The usefulness of XCDCC is varied. Near 
future apphcations of XCDCC include ^^Be(p,p')^^Be*, 
^^C(^"^Be,^''Be-|-n), and the various modes of ^B breakup. 
A comparison with breakup data to specific core states 
will allow accurate checks of the validity of the structure 
model for the projectile and thereby improve the compar- 



ison between structure calculations and experiment. We 
will also be able to look at detailed breakup spectra and 
analyze resonances built on excited core states, which 
could not be studied using previous reaction models. 

The generality of the formalism for including core de- 
grees of freedom in XCDCC shows that one can go be- 
yond the simple particle-rotor model for the projectile. 
Given the radial overlap wave functions between the pro- 
jectile and the core states, and the core elastic and in- 
elastic transitions for the core-target scattering, a more 
microscopic picture of the reaction can be obtained, in- 
cluding important dynamical effects which can have a 
significant effect on the reaction observables. 
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APPENDIX A: DETAILS ON THE THREE-BODY 
COUPLING POTENTIALS 

Here we give the details of the derivation for Eq. I|23|l . 
The expansion in Q, Eq. (|18|l . and A, Eq. (|20|l . gives us 
three spherical harmonics which depend on the R, f and 
^ degrees of freedom. The expansion in K , Eq. (|21|l . gives 
us two spherical harmonics in R and r, which combine 
with the spherical harmonics from the Q expansion to 
give two new multipole orders, A and A', which relate 
to the projectile-target and core-valence relative motion 
respectively: 



Ciffc(^)CA;.(A) -^Ca^(H)A2(-1) 





A A \ 




A 

















(Al) 



(r)-^CA'.'(r)A''(-ir' 



K Q-X A' 

— k q — fi —UJ' 



K Q-X A' 




(A2) 



We then combine the three spherical harmonics, for the 
three coordinates, to form a tripolar spherical harmonic 
with zero total angular momentum, which is a tensor 
of rank zero, since the matrix elements of the potential 
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operator must conserve total angular momentum: 

(_1)A'-Q /A A' Q 



A 



y.{CA{R)-{CA,{f)®CQ{0}K) ■ 

(A3) 

Using these definitions, one can write the potential oper- 



ator as Eq. (EH). 



In order to calculate the matrix elements, we use the 
projectile coupling order defined in Eqs. ©fe©- The 
matrix elements of the tripolar spherical harmonic of 
rank zero is: 



J 



m{ls)j, I; Jp])J, Jt; Jt\{Ca{R) ■ {CA,{f) ® Cq{0}a)\{L'[{1's)j' ,1'; J'p])J', J[- Jj 



A L L' 




A' / V 




Jp J'p A 
L' L J 



J f A'l 
I' I s 



' Jp J'p A 
3 f A'}{I\\Cq{0\\I') 
I I' Q 



r 



(A4) 



By rearranging the summations, one can write the tran- 
sition potentials as defined in Eqs. (|^ - l?7| . 

The reduced matrix element for the core degree of free- 
dom is defined using the Bohr & Mottelson definition: 



The radial integrals are 







{I II Cq{0 II /') = I{I\CQiO\I')/iI'K'^Qq\IK,) . (A5) and the couplings are 



(B3) 



We have left this matrix element separate since the 
derivation of the coupling potential is independent of the 
choice for the core inelastic excitation mode. The specific 
form for the matrix element within a rotational model is 
presented in Appendix C. 



APPENDIX B: Q = LIMIT 

If one sets Q = in the transition potential defined in 
Eqs. (PlH?7jl . this corresponds to the limit that the core 
cannot dynamically excite during the reaction, hence the 
core spin in the initial state is the same as the core spin in 
the final state. In the Q = limit, the coupling potentials 
are the same as for single particle breakup as defined in 
[58l |. since A = K, 

u'JAR) =LL'jpj'p{~iY-+'' Y.^-ifk^ 

K 





L 




{Jp 


J'p 




(o 







u 


L 





(Bl) 

where the form-factors are the sum over the coupled 
channels states of the radial integrals and couplings. 



p^T'' =Sii' i-iy' Kjj'ii' 

K I l'\ (Jp J'p K\ (j f K 

on / J i\\i' I s 



X 



(B4) 



One can see from this limit that once the form-factors 
are calculated, the coupling potentials have same form 
as with core excitation. We associate this limit with a 
reaction model where the core cannot dynamically ex- 
cite due to the interaction with the target. This limit 
may still include couplings between the core and valence 
particles which can excite/de-excite the core. 



APPENDIX C: ROTATIONAL MODEL FOR THE 
CORE 

In the formalism for XCDCC, we did not assume a 
model for the core- valence coupling. Here we evaluate the 
XCDCC equations using a rotational model for the core. 
For the radial dependence, we use a deformed Woods- 
Saxon potential 



Kc(r,0 = 



V 



(CI) 



1 + exp [(r - Rm/av] ' 
where the radius depends on the quadrupole deformation: 



i?(0=i?l/(l + /32r2o(0) 



(C2) 
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The core- valence potential are numerically calculated by 
expanding in multipoles, as defined in |58l |. 



(C3) 



In the rotational model, the eigenstates of the core, ipi{£,), 
are proportional to the rotational matrices, 



(C4) 



The matrix element, for the core degrees of freedom, in 
the coupling potentials is then 



{iWCgiO II /') = (-1/'+^^//' 



Q I r 

K, 



(C5) 



where Kc is the projection of the core spin, /, on the 
body-fixed z-axis. 



[1 

[2: 

[3] 
[4] 
[5] 
[6] 
[7] 
[8] 

[9: 
[lo: 

[11 

[12 

[13 
[14 
[15 
[16 
[17 

[18 

[19 

[20 

[21 

[22' 
[2i 

[24 

[25 
[26 
[27 
[28; 

[29 

[30 

[31 

[32; 



B. Jonson, Phys. Rep. 389, 1 (2004). 
J. S. Al-Khalili and F. M. Nunes, J. Phys. G 29, R89 
(2003). [33 
T. Nakamura et al, Phys. Rev. Lett. 83, 001112 (1999). 
J. H. Kelley et al, Phys. Rev. Lett. 74, 30 (1995). 
N. A. Orr et al, Phys. Rev. C 51, 3116 (1995). [34 
V. Maddalena et al, Phys. Rev. C 63, 024613 (2001). 
T. Aumann et al, Phys. Rev. Lett. 84, 35 (2000). [35 
V. Maddalena and R. Shyam, Phys. Rev. C 63, 
051601(R) (2001). 
N. Fukuda et al, Phys. Rev. C 70, 054606 (2004). [36 
B. M. Young et al, Phys. Rev. C 49, 279 (1994). [37 
K. Alder and A. Winther, Electromagnetic Excitation 
(North-HoUand, Amsterdam, 1975). [38 

G. Baur, K. Hencken, and D. Trautmann, Prog. Part. 

Nucl. Phys. 51, 487 (2003). [39 
T. Motobayashi et al, Phys. Lett. B264, 259 (1991). 
T. Kikuchi et al, Phys. Lett. B391, 261 (1997). [40 

B. Davids et al, Phys. Rev. Lett. 81, 2209 (1998). 

F. Fleurot et al, Phys. Lett. B615, 167 (2005). [41 

H. Esbensen and G. F. Bertsch, Nucl. Phys. A600, 37 
(1996). ' [42 

C. H. Dasso, S. M. Lenzi, and A. Vitturi, Nucl. Phys. 
A611, 124 (1996); A639, 635 (1998). 

F. M. Nunes and L J. Thompson, Phys. Rev. C 57, R2818 [43 

(1998) . 

H. Esbensen and G. F. Bertsch, Phys. Rev. G 59, 3240 [44 

(1999) . 

F. M. Nunes and L J. Thompson, Phys. Rev. G 59, 2652 [45 
(1999). [46 
P. Gapel and D. Baye, Phys. Rev. G 71, 044609 (2005). 
K. Ogata, S. Hashimoto, Y. Iseri, M. Kamimura, and M. [47 
Yahiro, arvix.org/nucl-th/0505007 (2005), submitted to 
Phys. Rev. G. [48 
H. Esbensen, G. F. Bertsch, and K. A. Snover, Phys. [49 
Rev. Lett. 94, 042502 (2005). [50 
J. J. Kolata et al, Phys. Rev. G 63, 024616 (2001). 

D. Gortina-Gil et al, Nucl. Phys. A720, 69 (2003). 

F. M. Marques et al, Phys. Lett. B476, 219 (2000). [51 

J. A. Tostevin, S. Rugmai and R. G. Johnson, Phys. Rev. 

G 57, 3225 (1998). [52 

E. Garrido, D. V. Fedorov and A. S. Jensen, Phys. Rev. 

G 59, 1272 (1999). [53 
Y. Sakuragi, M. Yahiro, and M. Kamimura, Prog. Theor. 
Phys. Suppl. 89, 136 (1986). [54 
J. A. Tostevin, F. M. Nunes, and L J. Thompson, Phys. 
Rev. G 63, 024617 (2001). [55 
T. Matsumoto, E. Hiyama, K. Ogata, Y. Iseri, M. 



Kamimura, S. Ghiba, and M. Yahiro, Phys. Rev. G 70, 
061601 (R) (2004). 

T. Egami, K. Ogata, T. Matsumoto, Y. Iseri, M. 
Kamimura, and M. Yahiro, Phys. Rev. G 70, 047604 

(2004) . 

M. Kamimura, T. Matsumoto, E. Hiyama, K. Ogata, Y. 
Iseri, and M. Yahiro, AIP Gonf. Proc. 791, 174 (2005). 
T. Matsumoto, T. Egami, K. Ogata, Y. Iseri, M. 
Kamimura, and M. Yahiro, arvix.org/nucl-th/0511059 

(2005) , submitted to Phys. Rev. G. 

G. A. Bertulani, Phys. Rev. Lett. 94, 072701 (2005). 

E. O. Alt, B. F. Irgaziev, and A. M. Mukhamedzhanov, 
Phys. Rev. G 71, 024605 (2005). 

K. Rusek and K. W. Kemper, Phys. Rev. G 61, 034608 

(2000) . 

M. Takashina, S. Takagi, Y. Sakuragi, and Y. Iseri, Phys. 
Rev. G 67, 037601 (2003). 

A. M. Moro, R. Grespo, F. M. Nunes, and I. J. Thomp- 
son, Phys. Rev. G 67, 047602 (2003). 

N. G. Summers and F. M. Nunes, Phys. Rev. G 70, 
011602(R) (2004). 

B. Davids, Sam M. Austin, D. Bazin, H. Esbensen, B. M. 
Sherrill, I. J. Thompson, and J. A. Tostevin, Phys. Rev. 
G 63, 065806 (2001). 

J. Mortimer, I. J. Thompson, and J. A. Tostevin, Phys. 
Rev. G 65, 064619 (2002). 

N. G. Summers and F. M. Nunes, J. Phys. G 31, 1437 
(2005). 

A. Shrivastava et al, Phys. Lett. B596, 54 (2004). 

N. G. Summers, F. M. Nunes, and I. J. Thompson, Phys. 

Rev. G (R) (to be published). 

P. G. Hansen and J. A. Tostevin, Annu. Rev. Nucl. Part. 
Sci. 53, 219 (2003). 

J. A. Tostevin, J. Phys. G 23, 735 (1999). 
J. A. Tostevin, Nucl. Phys. A682, 320c (2001). 
J. A. Tostevin, D. Bazin, B. A. Brown, T. Glasmacher, 
P. G. Hansen, V. Maddalena, A. Navin, and B. M. Sher- 
rill, Phys. Rev. G 66, 024607 (2002). 
R. Shyam and P. Danielewicz, Phys. Rev. G 63, 054608 

(2001) . 

F. M. Nunes, I. J. Thompson, and R. G. Johnson, Nucl. 
Phys. A596, 171 (1996). 

F. M. Nunes, J. A. Ghristley, I. J. Thompson, R. G. John- 
son, and V. D. Efros, Nucl. Phys. A609, 43 (1996). 
D. Ridika, M. H. Smedberg, J. S. Vaagen, and M. V. 
Zhukov, Nucl. Phys. A628, 363 (1998). 
T. Tarutina and M. S. Hussein, Phys. Rev. G 70, 034603 
(2004). 



14 



[56] P. Batham, I. J. Thompson, and J. A. Tostevin, Phys. 

Rev. C 71, 064608 (2005). 
[57] J. A. Christley and I. J. Thompson, Comp. Phys. Comm. 

79, 143 (1994). 
[58] I. J. Thompson, Comput. Phys. Rep. 7, 3 (1988). 
[59] M. Moshinsky, Nucl. Phys. 13, 104 (1959). 



[60] A. Ozawa ei al. Nucl. Phys. A693, 32 (2001). 
[61] M. G. Saiut-Laureut ei al, Z. Phys. A 332, 457 (1989). 
[62] J. S. Al-Khahh, J. A. Tostevin, and I. J. Thompson, 
Phys. Rev. C 54, 1843 (1996). 



